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ABSTRACT 

We study whether integrated optical spectroscopy of a disk galaxy can be used to infer the mean, or 
characteristic gas-phase oxygen abundance in the presence of systematic effects such as spatial abun- 
dance variations, contributions to the integrated emission-line spectrum from diffuse-ionized gas, and 
dust attenuation. Our sample consists of 14 nearby disk galaxies with integrated spectrophotometry, 
and observations of more than 250 individual H ii regions culled from the literature. We consider both 
theoretical and empirical strong- line abundance calibrations based on the R23 = ([O 11] + [O iii])/H/3 
parameter. We find that the integrated oxygen abundance correlates well with the gas-phase abun- 
dance measured at a fixed galactocentric radius, as determined by the H ii-region abundance gradi- 
ent. The typical scatter in the correlation is ±0.1 dex, independent of the abundance calibration, or 
whether the observed integrated emission-line fluxes, the reddening-corrected fluxes, or the emission- 
line equivalent widths are used. Integrated abundances based on the observed fluxes or equivalent 
widths, however, are susceptible to additional systematic effects of order 0.05 — 0.1 dex, at least for 
the range of reddenings and stellar populations spanned by our sample. Unlike the integrated R23 
parameter, we find that the integrated [N 11] /Ha and [S 11] /Ha ratios are enhanced with respect to 
line-ratios typical of H 11 regions, consistent with a modest contribution from diffuse-ionized gas emis- 
sion. We conclude that the R23 parameter can be used to reliably measure the gas-phase abundances 
of distant star-forming galaxies. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: spiral 



1. INTRODUCTION 

Improving constraints on the chemical evolution 
histories of star-forming galaxies requires abundance 
measurements with well-understood systematic uncer- 
tainties. Quantifying these uncertainties is particu- 
larly important in light of recent results which show 
that typical star-forming galaxies undergo a rela- 
tively modest amount of chemical enrichment since 
z ^ 1 (J . Mou s takas et al. 2006, in prepa ration: 
i Lillv ^t alJ 120031: iMaier et alJ r2 004: Kobulnick v et ijil.l 
12003: Kobulni ckv fc Kewlevl '2004: S avadio et alJ i2005). 

In the nearby universe, galaxies are resolved into in- 
dividual H II regions, enabling a direct measurement 
of the gas-phase abundance at discreet spatial posi- 
tions. By comparison, spectra of distant galaxies are 
a luminosity-weighted average of star-forming regions 
spanning a range of physical conditions, stars, and diffuse 
gas. In this paper we investigate whether spatially unre- 
solved spectroscopy of spiral galaxies encompassing all, 
or a significant fraction of the light, reliably measures the 
gas-phase abundance as determined from spectroscopy of 
individual H 11 regions. 

Our analysis focuses on disk galaxies for two reasons. 
First, luminous star- forming disk galaxies at high red- 
shift are more likely to be observed compared to low- 
luminosity dwarf galaxies. Second, several physical prop- 
erties of disk galaxies may systematically bias abundance 
measurements based on integrated spectroscopy. For ex- 
ample, whereas spiral galaxies exhibit radial abundance 
variations, to first order, dwarf galaxies are chemically 
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homogenous re.g.. lKobulnickv SkillriT^rmfl^) . There- 
fore, we anticipate the largest discrepancy between inte- 
grated and spatially resolved abundances to arise in a 
sample of disk galaxies. 

One of the most striking and characteristic proper- 
ties of disk galaxies is the tendency for their centers to 
be more metal rich than their outskirts. Although the 
exact physical processes that give rise to radial abun- 
dance gradients are still a matter of debate l)Molla et alJ 
119961 and references therein), it has been shown that 
the vast majority of near by d isk galaxies, if not all, pos- 
sess radial gradients (e.g..lVil a-Costas & Edmunds 199^ 
iZaritskv et alJll994t Ivan Zee et al.lll998^1 . Since distant 
galaxies may have steeper abundance gradi ents, indica- 
tive of their earlier evolutionary state fe.g.. iMolla et alJ 
it is important to quantify how variations in the 
slope of the abundance gradient affects the abundances 
we infer from integrated spectroscopy. 

Another common, if not ubiquitous physical property 
of spiral galaxies, including the Milky Way, is the diffuse, 
low surface-brightness phase of the interstellar medium 
called the diffuse- io nized gas (DIG), or the diffuse- 
ionized medium (e.g., Rcvnolds 1990: Haffncr ct al. ll999t 
[Ferguso n ct al. 199 6; Wang et al. 1997; Grcenawal t et alJ 
1998; IZurita et alJ i2000: Thilkcr ct al. 2002). Of or- 
der 30 — 50% of the integrated Ha luminosity in spi- 
ral galaxies can be attributed to the DIG, making it 
an energetically important component of the interstellar 
medium of disk galaxies (Calzetti et al. 2004, and ref- 
erences therein). Since low- ionization line-ratios such as 
[N 11] /Ha and [S 11] /Ha in the DIG are enhanced relative 
to ra tios involving high-excitation lines such as [O iii]/H/3 
(e.g.. lMa,rtiu1IT99lmm^es Wa ,lterbos"200,?). we want 
to determine what effect, if any, the integrated DIG emis- 
sion of a disk galaxy has on the abundance we measure 
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from its integrated spectrum. 

Finally, it is important to explore how variations in 
dust attenuation might bias integrated abundance mea- 
surements. Since the Ha line becomes inaccessible to 
ground-based optical spectroscopy above z « 0.4, it may 
not always be possibl e to correct the nebular lines for 
dust reddening (e.g. 



„ Kobulnickv et al 

e.g., iLiang et alJ 12004 and iMaier et al 
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but see. 
There- 
fore, we compare the abundances derived from the ob- 
served line fluxes, the reddening-corrected fluxes, and the 
emission-line equivalent widths, which, to firs t order, are 
insen sitive to the effects of dust ()Kobulnickv fc PhiUips 
I200I . In the presence of all these systematic effects, we 
want to determine whether we can measure the "inte- 
grated" gas-phase abundance of a disk galaxy, and, if so, 
how we can physically interpret that measurement. 

Previous efforts investigating the integrated abun- 
dances of spiral galaxies found that integrated spec- 
troscopy provides a fairly reliable measure of the gas- 
phase abundance ijKobulnickv et al.lI999t IPilvugin et al.l 
12004a,) . However, these studies synthesized integrated 
spectra by co-adding observations of individual H 11 
regions; therefore, they could not investigate the sys- 
tematic effects of DIG emission and dust reddening. 
Recently, we have obtained spatially unbiased (inte- 
grated) spectroscopy for severa l hundred nearby galaxies 
l|Moustakas fc Kennicutt 120061 hereafter MK06). Here, 
we combine these data with spectroscopy of individ- 
ual H II regions culled from the literature to investi- 
gate whether chemical analysis methods can be used to 
estimate the metallicities of spatially unresolved star- 
forming galaxies. 

2. DATA 

2.1. Sample Selection and Emission- Line Measurements 

Investigating the integrated nebular abundances of 
disk galaxies quantitatively poses several observational 
challenges. First, because of practical considerations, 
H ii-region abundance measurements are generally only 
available for nearby galaxies subtending relatively large 
(> 2') angular diameters on the sky. At the same time, 
obtaining an integrated spectrum of an object that is 
larger than the typical length of a long-slit (3' — 5') can be 
very difficult. Nevertheless, our survey (MK06) includes 
among the first integrated spectra of spiral galaxies with 
measurable radial abundance gradients. 

As part of a larger effort to study extra-galactic star 
formation, we have obtained integrated, optical (3600 — 
6900 A) spectrophotometry at ~ 8 A FWHM resolu- 
tion for a diverse sample of 417 nearby (< 150 Mpc) 
galaxies (MK06). Our survey implements the lKennicutti 
p992i|1 long-slit drift-scanning technique to obtain inte- 
grated spectra of large galaxies at the spect ral resolution 
afforded by a 2'.'5 x 200" long-slit (see also l.Tansen et al.1 
[2000; Gavazzi et al. 2004). The size of each galaxy at 
the B25 mag arcsec"^ isophotc dictates the parameters 
of the drift-scan, and the width of the extraction aper- 
ture. The resulting spectra typically include 80 — 100% 
of the optical light. These observations are ideally suited 
to our analysis because they are analogous to obtaining 
a normal (spatially fixed) optical spectrum of a distant 
galaxy, which subtends a corresponding smaller angular 
diameter on the sky. 



In MK06 we present the fluxes and equivalent widths 
(EWs) of the strong nebular emission lines, includ- 
ing [O 11] A3727, [O III] A5007, [N 11] A6584, Ha, 
and H/3, for the complete sample of galaxies. Here- 
after, we assume that [O ill] A4959 is 0.34 times 
the [O I II] A5007 intensity, as determined by atomic 
physics ijStorev fc ZeippenI 12000(1 . Unfortunately, the 
temperature-sensitive [O ill] A4363 line is typically not 
detected in our spectra, which is not surprising given that 
this line becomes vanishingly weak in metal-rich galaxies 
(e.g., Kennicutt et al. 2003 ); therefore, we must rely on 
so-called strong-line calibrations to estimate the oxygen 
abundance (see ^3.1|l . 

The nebular line measurements presented in MK06 in- 
clude a careful treatment of underlying stellar absorp- 
tion, which is among the most important systematic ef- 
fects to consider in any emission-line abundance analysis. 
In spectroscopic studies of individual H 11 regions, the 
Balmer lines are typically corrected for a constant ~ 2 A 
of stellar absorption, corresponding to the amount of ab- 
sorption exp ected for a yo ung (< 20 Myr) stellar pop- 
ulation ('e.g.. lMcCall et alj fl985: Olofsson 1995). How- 
ever, in general, the integrated, luminosity-weighted stel- 
lar population of a galaxy will differ considerably from 
that of a single star-forming region. In particular, the 
assumption of a constant amount of stellar absorption 
for Ha and H/3 may break down, since different stel- 
lar populations will dominate the integrated light at red 
(~ 6500 A) and blue (~ 4800 A) wavelengths, according 
to the particular star-formation history of each galaxy. 
Moreover, the exact correction depends on spectral reso- 
lution, since the nebular lines are typically narrower than 
the underlying absorption. 

To address these issues, in MK06 we use population 
synthesis to model the integrated absorption-line spec- 
trum of each galaxy in our sample. We fit a non- 
negative linear combinatio n of seven instantaneous-burst 
iBruzual fc Chariot 1)20031) population synthesis models, 
convolved to the instrumental resolution, with ages rang- 
ing from 5 Myr to 12 Gyr. Including a simple treatment 
of dust reddening, this method results in precise sub- 
traction of the stellar co ntinuum underlying the nebular 
emiss i on lines (see also |Ho et al.lll997t iTremonti et alJ 
120041: ICid Fernandes et all 120*05^ . For the subset of 
galaxies considered in this analysis (see below), the H/3 
absorption ranges from 2.8 to 4.9 A, with a mean value 
of 3.9 ± 0.5 A. The corresponding correction for Ha is 
1.8 - 3.0 A, with a mean of 2.5 ± 0.3 A. 

We define our sample according to two selection cri- 
teria. First, we restrict the objects observed by MK06 
to spiral galaxies with well-detected (30-) Ha and H/3 
emission lines. We also reject galaxies whose integrated 
emission-line spectrum is clearly dominated by an active 
nucleus (e.g., NGC 1068). Next, we conduct a search 
of the literature for spectroscopy of individual H 11 re- 
gions in the remaining objects."^ Our search results in 
data on 252 H 11 regions in 14 disk galaxies. The me- 
dian number of H II regions per galaxy is 15. One ob- 

^ We only consider observations that include the [O 11], [O III], 
and H/3 nebular lines, since these are the most likely lines to 
be measured in distant galaxies. This requirement eliminates 
the imaging spec trophotometr i c data of the galaxies in our sam- 
ple presented by iDutil fc Rovl {19991) . IMartin fc Rol 119921) . and 
llviartin fc BellevI 119971) . arnong others. 
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ject, NGC 4713, has fewer than six H ii regions, which 
IZaritskv et alJ ^^1994^ advocate as the minimum number 
of regions needed to robustly determine the abundance 
gradient (see alsofbutil & Rov 2001). Nevertheless, we 
retain NGC 4713 in our sample, bearing in mind that 
the derived abundance gradient may be subject to addi- 
tional systematic uncertainty. To place the size of our 
sample in context, these 14 objects account for ~ 1/4 of 
the total numb er of galaxies with mea surable abundance 
gradients fe.g.. lPilvugin et al]l2004bl) . 

In Table ^ we list the galaxies in our sample and 
the global properties relevant to our analysis. Column 
(2) gives the morphologica l type of each object from 
Ide Vaucouleurs et al.l (|1991, hereafter RC3). Column (3) 
lists the absolute B-band magnitude, based on the Galac- 
tic extinction- corrected i?-band mag nitude [Ry = 3.1; 
iSchlegel et al] ll998: O' DonnellHigpT and distance tabu- 
lated in MK06. Column (4) lists P25, the radius of the 
semi-major axis at the B25 mag arcsec"^ isophote from 
the RC3. In columns (5) and (6), respectively, we list the 
photometric inclination a ngle and posit ion angle, as de- 
termined in the X^-band ijJarrett et al.l l2000. 2003). To 
compute the inclination angle we use the observed major- 
to-minor axis ratio at th e Ks,20 mag arcsec"^ isophote 
ijJarrett et alJl2000l 1200 3) , and assume that the axial ra- 
tio of a system viewed edgc-on is 0.2 (TuUy ct al. 1998). 
For NGC 5194, however, we us e the kinematic inclina- 
tion angle determined by iTuUvl l)1974) . Finally, in col- 
umn (7) we give the Ha surface brightness, I](Ha) = 
4.5 X 10^'*/(IIa)/p|5 in erg s~^ pc~^, where /(Ha) is the 
Ha flux in erg s~ cm~^ from MK06, corrected for ex- 
tinction as described in W6.1[ The remaining columns in 
Tableware described in H3.ll 

For each H 11 region drawn from the literature, we 
tabulate the coordinates relative to the galactic nu- 
cleus, and the stellar absorption and reddening-corrected 
[O 11] A3727, [O III] A5007, and H/3 emission-line flux 
measurements and uncertainties. As before, we set 
[O III] A4959 = 0.34 x [O iii] A5007. Using the optical 
disk radius and galactic position and inclination angles 
given in Tabled we also compute the de-projected galac- 
tocentric radius of each H 11 region, normalized to p25. In 
ti3.2l we use these data to determine the radial abundance 
gradients for the galaxies in our sample. 

2.2. Sample Properties and Selection Biases 

In this section we discuss some of the sample selec- 
tion biases. As alluded to in H2.ll the size of the sam- 
ple is determined by simultaneously requiring integrated 
spectroscopy and spectroscopy for a sufficient number of 
individual H 11 regions. In terms of i?-band luminos- 
ity, the sample is generally biased toward bright galax- 
ies, since lower-luminosity galaxies tend not to exhibit 
abundance gradients. We note, however, that our sample 
spans more than an order-of-magnitude in _B-band lumi- 
nosity, from Mb — —18.2 to —21.0 mag. For comparison, 
the absolute B magnit udes of the disk galaxies studied by 
IPilvugin et al.l l|20041j) , which is the most comprehensive 
compilation of galaxies with measured abundance gradi- 
ents to date, range from —17.8 to —21.6 mag. The distri- 
bution of inclination angles in our sample is fairly broad, 
ranging from 11° in NGC 3344, to 72° in NGC 3198; 
the mean inclination angle is 41° ± 19°. In 2|we test 
whether variations in inclination along the line-of-sight 



systematically biases the integrated abundances. 

The distribution of morphological types in our sam- 
ple is weighted toward late- type galaxies: 85% of the 
sample is type Sbc or later, with only two early-type 
objects, NGC 3351 (Sb) and NGC 4736 (Sab). This 
bias arises because of observational and physical se- 
lection effects which tend to prefer late-type galaxies. 
For example, the number of early-type galaxies with 
measured abundance gradients is relatively limited, be- 
cause H II regions in early-type galaxies are rarer and 
have lower luminosities and surface brightnesses than 
corresponding H 11 regions in late-type galaxies (e.g., 
IKennicuttil988|) . Previous studies have found that early- 
type galaxies are, on average, more metal rich, and have 
shallo wer a bundance gradients than late-type galaxies 
(e.g.. iGaT^ett fc Shields 1987; Gey fc Kennicutt 1993 
IZaritskv et^ ll994: Dutil & Rov) 119991) . We note, how- 
ever, that our sample includes galaxies that are as lu- 
minous and metal-rich as ot her well-known early-typ e 
disk galaxies such as M 81 l|Garnett fc Shiclds"'l987'). 
A related selection effect is that in the integrated spec- 
tra of early-type galaxies the equivalent widths of the 
nebular lines, if present, tend to be signiflcantly lower 
than the corresponding emission-line equivalent widths 
in late- type galaxies. This observation can be attributed 
to the lower star-formation rate p er unit luminosity or 
mass in early-type ga laxies (e.g.. ICaldwell et all 119911: 
iKennicutt et al.lll994|) . 

Despite the relatively small number of early-type 
galaxies, we argue that our sample is still representa- 
tive for the kind of analysis we are conducting, and that 
our conclusions are not strongly affected by this mor- 
phological bias. In particular, the objects most likely 
to be studied at high redshift should have higher gas 
resevoirs, high specific star- format ion rates, and large 
emission-line equivalent widths. For example, in an 
analysis of 64 sta r-forming galaxies at 0.3 < z < 0.8, 
iKobulnickv eTal] l|2003^) found an average H/3 em ission- 
hne EW of 21 ± 14 A, while iPettini et all 1)200 ID mea- 
sured H/3 EWs in excess of ~ 20 A in a sample of five 
Lyman-break galaxies at z ~ 3. Finally, since early-type 
galaxies tend to have shallow abundance gradients, small 
amounts of dust, and a weak or absent DIG phase, we 
do not expect the integrated abundances of early-type 
galaxies to differ significantly from the abundances in- 
fered from spectroscopy of individual H 11 regions. 

3. ANALYSIS 

3.1. Integrated Abundances 

We use the emission-line measurements from MK06 
to compute the integrated oxygen abundances of the 
galaxies selected in ii2.1l A principal assumption of 
this method is that the observed emission lines arise 
via photoionization from massive stars. Consequently, it 
is important to assess any possible contamination from 
nuclear activity in our integrated emission-line spectra. 
iHo et all l)1997t) have obtained nuclear spectra of all the 
galaxies in our sample; they find that 50% (7/14) have 
H ii-hke nuclei, 36% (5/14) are LINERs or transition 
objects, and 14% (2/14) exhibit Seyfert 2 nuclear line- 
ratios. In all cases, however, we find that the nucleus con- 
tributes < 3% to the integrated Ha flux. Therefore, we 
can apply standard techniques to derive chemical abun- 
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dances from the observed emission-line spectra. 

From the observed Ha/H/3 ratio, we estimate the neb- 
ular reddening assuming an intrinsic ratio of 2.86 and the 
IP'Donncll (1994) Milky Way extinction curve (listed in 
col. [8] of Tabic nj. The distribution of reddening val- 
ues for the sample ranges from 0.08 to 0. 69 mag, with a 
mean value of 0.31 ± 0.17 mag. Mousta kas et al.l l|2006n 
show that adopting a simple Milky Way extinction law 
works remarkably well at accounting for the effects of 
dust on the nebular lines from integrated optical spectra 
(see also Kewl ev ct al. 2002), despite the wide variation 
in attenuation l aws pre dicted by theoretical mo dels (e.g., 
IWitt fc Cordoni.200nt.Fischera Donitall200l. 

To investigate the systematic effect of dust atten- 
uation, we inter-compare the abundances determined 
using the observed (un-corrected) emission-line fluxes, 
the re ddening-corrected flux e s, an d the emission-line 
EWs. iKobulnickv fc Phillip^ l)2003j) have suggested re- 
placing fluxes with EWs when determining the integrated 
emission-line abundances of galaxies, since, to first or- 
der, EWs are insensitive to dust attenuation. However, 
the EW method for deriving abundances is sensitive to 
intrinsic variations in the stellar populations, and to dif- 
ferential extinction between the lines and the continuum, 
as we discuss in 21 

We explore the effect of choosing a particular 
strong-line c alibration by com puting abundances us- 
ing both the iMcGaughl l|l99lL hereafter M91) and the 
iPilvugin fc ThuanI l|200a hereafter PT05) calibrations of 
the i?23 = ([O II] A3727 -h [O lii] AA4959, 5007)/H/3 pa- 
rameter. We focus on the R23 parameter as an abun- 
dance diagnostic because the oxygen and H/3 lines are 
among the strongest optical emission lines, and they 
are observable in ground-based optical spectra to z 
1. Both calibrations account for variations in ioniza- 
tion parameter at fixed metallicity, either through the 
[O iii]/[0 II] ratio (M91), or the [O iii]/([0 ii]-|-[0 ill]) 
ratio (PT05). The M91 calibration, as parameterized by 
IKobulnickv et all (^999), is based on an extensive grid of 
photoionization models spanning a wide range of metal- 
licity and ionization parameter. By comparison, PT05 
have constructed a purely empirical calibration, based 
on a carefully selected sample of H 11 regions with well- 
measured electron temperatures. Finally, we note that 
although i?23 is double- valued as a function of abundance 
(see, e.g., M91), all the galaxies in our sample lie on the 
upper branch of the relation between R23 and (0/H). 

On average, the M91 abundance scale is shifted to 
higher values relative to the PT05 scale by ~ 0.5 dex, 
or a factor of 3 (see Fig. EJ. In 21 we show that se- 
lecting among the M91 and PT05 calibrations has a sig- 
nificant effect on the absolute oxygen abundances, and 
a modest effect on the slope of the derived abundance 
gradients, but does not significantly alter our conclu- 
sions. Using a diffe rent th eoretical calibration in place of 
M91 (e.g. , Zaritsk v'^ral]ll9 94: Kcwlcv & Dooita 2002; 
IKobulnickv & Kcw levl200ir has a negligible effect on our 
conclusions. 
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Fig. 1. — Oxygen abundance vs. de-projected galactocentric ra- 
dius, p, normalized to p25, the disk radius at the B25 mag arcsec"^ 
isophote. Oxygen abundances for each H II region have been com- 
puted using both the McGaugh (1991, M91) (circles) and the 
IPilvugin fc ThuanI 12005 . PT05) (squares) abundance calibrations. 
The solid line is the best fit through the M91 abundance distribu- 
tion, and the dashed line is the corresponding fit through the PT05 
abundances. H II regions that have been excluded from the fits are 
shown as open symbols, with a vert ical line connecting the upper- 
and lower-branch abundances (see il3.2l for details). 

In Columns (9)-(ll) of Tablc^we list the three inte- 
grated abundances and statistical uncertainties derived 
using the M91 abundance calibration; columns (12)- 
(14) give the corresponding abundances based on the 
PT05 calibration. In every object except NGC 3351 and 
NGC 4321, the two dustiest galaxies in the sample, we 
detect the [O 11] , [O ill] , and H/3 emission lines with > 3(j 
significance. In NGC 3351 and NGC 4321, however, we 
derive an upper limit to the [O il] flux and EW, which 
translates into a lower limit on the oxygen abundance. 

3.2. H ii-Region Abundance Gradients 

We compute the radial abundance gradient of each 
galaxy in our sample using an unweighted linear least- 
squares fit to both the M91 and the PT05 H ii-region 
abundances. To quantify the uncertainty in the gradi- 
ent, we re-compute the fit 500 times, modulating the 
metallicity of each H 11 region by a Gaussian distribution 
with a width equal to the statistical uncertainty. TableEl 
and Figure ^ present the results of our fits. Columns 
(2)-(4) of Table (21 are based on the M91 abundance cal- 
ibration and list, respectively, the extrapolated central 
abundance at galactocentric radius p = 0, the charac- 
teristic abu ndance, defined as the oxygen abun dance at 
p = 0.4 P25 l|Zaritskv et al.llT99llGarnettll200l . and the 
slope of the abundance gradient in dex P2^. Columns 
(5)-(7) give the corresponding quantities derived using 
the PT05 abundance calibration. Finally, columns (8) 
and (9) give the number of H 11 regions used in the fit. 
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and the corresponding references to the literature, re- 
spectively. 

In Figure Q we plot the normalized galactocentric ra- 
dius, p/p25, versus oxygen abundance for the H ii regions 
in each galaxy. The circles and squares distinguish the 
M91 and PT05 abundances, respectively. We assume 
that all the H ii regions lie on the upper branch of the 
appropriate i?23-(0/II) relation, except as noted below. 
The solid line is the best fit through the M91 abundance 
distribution, and the dashed line is the corresponding fit 
through the PT05 H ii-region abundances. 

We exclude a handful of H II regions from the fits, 
shown in Figure ^ as open symbols without error bars, 
with a solid vertical line connecting the upper- and lower- 
branch i?23 abundances. Specifically, in NGC 1058 we 
exclude the t wo outermo st H ii regions, FGW1058G and 
FGW1058H ijFerguson e t al. 1998), because they are weU 
outside the main optical disk of the galaxy. We also 
excl ude the outermost H ii region in NGC 3198, -027- 
187 f Zaritskv et alj[l994|) . which has an extremely low 
ionization parameter, log([0 iii]/[0 ii]) = —0.76. Un- 
fortunately, the strong-line calibrations used here fail 
at reproducing the ionization properties of this object, 
which leads to the unphysical situation of having a lower- 
branch abundance that is several times higher than the 
corresponding upper-branch abundance. Finally, we also 
reject the anomolousl y low-metallicity H il region in 
NGC 3521, -033-118 ijZaritskv et aljll994j) . which has 
a negligible effect on the derived abundance gradient. 
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Fig. 2. — Slope of the radial abundance gradient determined 
from fitting to the [FilvuKin & Thuan 1 2005:i, PT05) abundance 
distribution vs. the corresponding slope derived from fitting to 
the McGauf^ IT99ll. M9f) H ll-region abundances (see Fig. [Tl . 
(Ins et) C orrelation between the characteristic oxygen abundances 
(see il3.2l and Table|2]for details). The dashed line shows the 0.5 dex 
mean offset between the M9I and PT05 abundance scales, while 
the solid line in both panels is the line-of-equality. 

In Figure El we compare the slope of the gradients 
derived using the M91 and PT05 abundance calibra- 
tions, and a comparison of the characteristic abundances 
(inset). The solid line in both panels is the line-of- 
equality, and the dashed line in the inset shows the me- 
dian 0.5 dex offset between the M91 and PT05 abun- 
dances. The two slopes generally correlate, although 
the PT05 gradients are systematically shallower than 
the corresponding gradients derived using the M91 oxy- 



gen abundances. The four most discrepant examples are 
NGC 3521, NGC 3344, NGC 1058, and NGC 4303. In 
all cases, the PT05 calibration predicts a much flatter 
relationship between galactocentric position and oxygen 
abundance, a result which has profound implications for 
our understanding of t he chemical evolution properties 
of sp iral galaxies (e.g.. lMolla et al.lir9 96: Pil vuein et alJ 
l2006j) . Exploring this result in more detail, however, is 
beyond the scope of this paper. 

4. RESULTS 

Figures 01 and 01 illustrate the principal result of this 
paper. Here we plot the characteristic oxygen abun- 
dance, as derived from the H ii-region abundance gra- 
dient (Fig. n Table O, versus our three estimates of 
the integrated abundance. We plot the integrated "ob- 
served", "reddening-corrected", and "EW" abundances 
using open circles, filled circles, and open squares, re- 
spectively, and the lower limits to the integrated abun- 
dances in NGC 3351 and NGC 4321 as arrows f^Oand 
Table [U. 

As an independent check on our results, in Fig- 
ure m we show the integra ted abundances determined 
by iKobulnickv et al.l l)1999()j against the charac teristic 
H ii-region abundances from'Zaritskv et ajj l)1994r [filled 
grey triangles). Kobulnickv et al. (1999) derive pseudo- 
integrated abundances for their sample of 22 galaxies, 
nine of which are common to our sample, by dividing 
each object into 5 — 13 equal-size radial zones and con- 
volving the observed Ha surface-brightness distribution 
with the reddening- c orrect ed H ii-region line-strengths 
from iZaritskv et all l)1994l) . We subtract 0.2 dex from 
these points to place them on the M91 abundance scale. 
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Fig. 3. — C haracteristic vs. integrated oxygen abundance, based 
on the IMcC^ ugh 1 1991) abundance calibration. Symbols differen- 
tiate between abundances computed using the observed emission- 
line fluxes (open circles), the reddening-corrected fluxes [filled cir- 
cles), and the emission-line equivalent widths (open squares). For 
comparison we plot the results from IKobulnickv et al. 11999) us- 
ing filled grey triangles, accounting for the 0.2 dex offset between 
the Zaritskv et al. 1 1994) and McGaush 1 1991) abundance scales. 
Lower limits to the oxygen abundance in NGC 3351 and NGC 4321 
are shown as arrows. The cross shows the median measurement 
uncertainty of the data. 

Examining Figures O and 01 we find that all three in- 
tegrated abundances correlate with the characteristic, or 
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spatially resolved abundance, independent of the abun- 
dance calibration used. A Spearman rank correlation test 
quantifies the significance of the observed correlations 
at > 99%. We conclude, therefore, that, to first order, 
the abundance inferred from the integrated emission-line 
spectrum of a galaxy is representative of the gas-phase 
oxygen abundance at p = 0.4 p25, even in the presence 
of an abundance gradient, contributions from DIG emis- 
sion, and variations in dust reddening. 
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Fig. 4. — Same as Fig. 13 but using the lPilvugin fc I'liuaJ ||200] 
abundance calibration. 



Table O presents the mean, standard deviation, and 
median of the residuals we measure between the in- 
tegrated and characteristic abundances. On average, 
the integrated abundance inferred from the observed 
(un-corrected) line-fluxes overestimates the characteris- 
tic abundance by 0.05 ± 0.06 dex or 0.09 ± 0.08 dex 
using the M91 or PT05 calibration, respectively. Ac- 
counting for dust reddening using a simple Milky Way 
extinction curve (O'Donnell 1994) removes the system- 
atic offset, but has almost no effect on the residual scat- 
ter. This result is consistent with other recent stud- 
ies, which find that the integrated emission lines of 
nearby star-forming galaxies appear to be attenuated by 
a simple foreground dust screen, and that the integrated 
B aimer decre ment is a robust indi c ator of the amount 
of dust fe.g.. ICalzetti et a,lJ ITflM I.Ta,nsen et abl l2noTt_ 



Kewlev et aLli2002t iDwita et alJl200a iMoustaLs et all 

The relatively small wavelength separation be- 
tween [O ii] A3727 and H/3 A4861, and the broad simi- 
larity of extinction and attenuation curves in the optical, 
additionally help to minimize the effects of dust. 

The integrated abundances determined from the 
emission-line EWs, by comparison, systematically under- 
estimate the characteristic abundance. The mean off- 
set is —0.06 ± 0.09 dex using the M91 calibration, or 
— 0.11±0.13 dex using the PT05 calibration. The median 
offsets are s lightly smaller, —0.04 dex and —0.09 dex, re- 
spectively. iKobulnickv fc Phillips! 1)20031) discuss in de- 
tail the systematic effects of using EWs to determine in- 
tegrated abundances. In particular, following their rec- 
ommendation, we have assumed that a = 1, where a is 
related to the ratio of the intrinsic continuum between 



3727 A and 4861 A, and the differential extinction be- 
tween the stellar continuurn and the emission lines (see 
IKobulnickv &: Phillips! 120031 their eq. [7]). Assuming 
a = 0.9, which is well wit hin the expected range of thi s 
parameter {a = 0.84 ± 0.3: "Kobulnick v fc Phillips!l200.'^ . 
increases the M91 integrated EW abundances by a me- 
dian -t-0.04 dex, and the corresponding PT05 abundances 
by H-0.06 dex. We conclude, therefore, that the inte- 
grated EW abundances are consistent with the character- 
istic gas-phase abundances, provided that an appropriate 
choice for a is made for the sample under consideration. 
We attribute the increased scatter in the integrated EW 
abundances relative to the abundances based on the line- 
fluxes to variations in a among star-forming galaxies. We 
note, however, that the measured scatter is comparable 
to the ±0.1 dex precision of the EW abundance technique 
iKobulnickv fc Phillips!l200l . 
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Next, we investigate the physical origin of the residual 
scatter between the integrated and characteristic abun- 
dances. In Figure |31 we plot the abundance residuals 
based on the M91 calibration, versus various physical 
properties of the sample: (a) the slope of the abun- 
dance gradient, (b) dust reddening, (c) inclination an- 
gle, and (rf) Ha surface brightness, E(IIq:), which char- 
acterizes the disk star-formation rate (SFR) per unit area 
ljKennicuttill998|) . For our sample I](IIa) spans a factor 
of ~ 15, comparable to the range in S(Hq;) cx E(SFR) 
found by^ennicutt ( 1998) for a larger sample of normal 
disk galaxies. The corresponding residual plots using the 
PT05 abundance calibration are shown in Figure El Ta- 
ble 01 presents the Spearman rank correlation coefRcients 
for each residual plot, and the probability that the vari- 
ables are uncorrelated. 

Examining Figures |S1 [HI and Table 0] we find that the 
residuals are not highly correlated with any of the vari- 
ables we have considered. We find a marginally signif- 
icant (> 80% probability) correlation between the ob- 
served and reddening-corrected M91 abundance resid- 
uals, and the galaxy inclination angle (Fig. |5j:;). We 
also find a marginally significant anti-correlation between 
these residuals and Tjijla) (Fig.|Sli). None of the residual 
correlations using the PT05 calibration are significant, 
although the increased scatter in the PT05 abundances 
(Fig-ffl Table|2Jl may be masking any underlying trends. 

The positive correlation between the abundance resid- 
uals and inclination angle, if it exists, is opposite of 
what we would have expected. The correlation indi- 
cates that the integrated abundance of an edge-on galaxy 
is higher than the corresponding metallicity of a less- 
inclined galaxy. Given the existence of radial abun- 
dance gradients (e.g.. Fig. P), we would have antici- 
pated the integrated spectrum of an inclined galaxy to 
be weighted more toward the outer, metal-poor H ii 
regions. This exact trend is shown by iTremonti et al.l 
p004 their Fig. 7) as a strong correlation between 
the residuals of the stellar-mass metallicity relation and 
inclination angle, based on an analysis of ~ 50,000 
star-forming galaxies in the Sloan Digital Sky Survey. 
ITremonti et al.l Ij2004f) find that edge-on galaxies may 
appear, on average, 0.2 dex (60%) more metal-rich 
than face-on galaxies. Since this systematic difference 
is comparable to the amount of chemical enrichment 



J. Mous- 



between the present day and z ~ 

takas et al. 2006, in preparati on; IKobulnickv et al.l 
120031 IKobulnickv fc Kewlevll2004j) . sample selection ef- 
fects that may be biased with respect to inclination an- 
gle should be carefully considered in look-back studies 
of galactic chemical e volution (e.g., TuUy-Fisher studies; 
iMouhcine et al. 2006). 

Since the strength of low-ionization line-ratios such 
as [N ii]/Ha and [S i i]/Ha correlate wit h the Ha sur- 
face brightness (e.g.. IWang et al.l 11998(1 . indicating a 
prominent DIG phase, the lack of a significant cor- 
relation between S(Ha) and the abundance residu- 
als in Figures and suggests that integrated 
abundance measurements are relatively unaffected by 
DIG emission. To explore this result in more de- 
tail, in Figure [7| we plot two emission-line diagnos- 
tic diagrams for our sample {squares), all the star- 
forming galaxies with integrated spectroscopy in the 
MK06 survey [diamonds), individual H ii regions drawn 
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Fig. 7. — Diagnostic diagrams showing the emission-line se- 
quences traced by our sample (squares), the full sample of star- 
forming galaxies from Moustakas & Kennicutt (2006) {diamonds), 
H II regions (grey points), and regions of diffuse-ionized gas in M 33, 
M 51a, and M 31 (crosses), as described in Jl] (a) Reddening- 
corrected [S ii] A6716/HO ratio vs. [N ll] A6584/Ha. (b) Reddening- 
corrected [O III] A5007/H/3 ratio vs. [O ii] A3727/H/3. 



from a dive rse sample o f disk and dwarf ga laxies 
" McGall et all IT985t IZaritskv et all IT991 



1 grey points; 



Ivan Zee et"aLlll998Hlzotov &: Thuanlll998D aildTeeions 
of DIG in three disk galaxies (crosses), NGC 0598=M 33, 
NGC 5194=M 51a, and NG C 0223=M 31 , base d on 
deep long-slit spectroscopy bv iGalarza et all l|1999() and 
iHoopes fc Walterbosi l|2003|) . 

Figure EJs shows that H ii regions and galaxies trace 
out similar emission-line sequences in the [N ii]/Ha 
versus [S ii]/Ha plane (corrected for reddening), while 
the DIG regions extend to much higher values of both 
[N ii]/Ha and [S ii]/Ha. IVIoreover, we observe that 
the integrated line-ratios preferentially lie at or above 
the uppe r envelope defined by the H ii regions (see also 
iLehncrt fc HeckmanI 11993) . The enhancement of low- 
ionization lines in the DIG is generally attributed to a low 
ionization parameter (e.g.. IPomgoreen fc Mathi^ 119941: 
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lMartirJfl997j) . Alternatively, recent self-consistent theo- 
retical models of aging star-forming regions reveal that 
old, low surface-brightness H il regions exhibit higher 
[N il]/Ha and [S ii]/Ha ratios relative to young, lumi- 
nous H II regions (M. Dopita 2006, private communi- 
cation). Regardless of its physical origin, our integrated 
spectra show a clear enhancement of low- ionization emis- 
sion relative to H ii regions. One implication of this 
result is that [N ii]- or [S ii] -based strong-line abun- 
dance diagnost ics applie d to integrated spectra of galax- 
ies (e.g., C hariot fc Longhctti 2001j Dcnicolo ct al. .200l 
iPettini fc P agcl 20^ may be moderately biased toward 
higher metallicitics. 

In Figure[7J) we plot the reddening-corrected [O iii]/H/3 
ratio versus [O ii] /H/3. In this diagram the emission- line 
sequences defined by galaxies, H ii regions, and regions of 
DIG generally overlap. Unhke [N Iij/Ha and [S ii]/Ha, 
lines of constant ionization parameter run nearly parallel 
to th e [O lii]/H/3 versus [O ii] /H/? emission-line sequence 
(e.s.. 'Ke wlev & DoDitall2002t) . We conclude, therefore, 
that measurements based on the globally averaged R23 
parameter give abundances that are fairly representative 
of the gas-phase oxygen abundances of galaxies. 

We conclude this section by testing one of the prin- 
cipal assumptions of the preceding analysis. Following 
IZarits kv et al. (199j), we have defined the characteristic 
abundance of a disk gala xy as the oxygen a b undance 
at p = 0.4 P25 (see also iGarnett et al.l Il997t iGarnetti 
l20f^ . However, given the form of the abundance gra- 
dient, we can derive the normalized galactocentric ra- 
dius, Pint/p25, where the integrated abundance equals 
the gas-phase abundance, assuming that the reddening- 
corrected integrated abundance is a perfect measure of 
the gas-phase abundance. For this calculation we only 
consider the M91 abundance calibration, which produces 
smaller uncertainties in the gradient slopes; we also ex- 
clude NGC 2541, which exhibits an increasing abundance 
with radius (although the slope is formally consistent 
with zero). 

We find that pint/p25 ranges from 0.2 in NGC 2903, 
to 0.8 in NGC 3344, with a mean value of 0.47 ± 0.15. 
Not surprisingly, Pint/p25 correlates with the slope of 
the abundance gradient, in the sense that the integrated 
abundances of galaxies with steeper abundance gradients 
are systematically higher than the integrated abundances 
of galaxies with shallower gradients. We conclude, there- 
fore, that the dominant source of scatter in the abun- 
dance residuals found in Figures EHHI (see also Table 
is due to assuming a fixed Pint/p25 = 0.4. Our analysis 
shows that while this is a reasonable average assump- 
tion, the exact radius does vary from object-to-object in 
a way that depends mildly on the slope of the abundance 
gradient. 

5. DISCUSSION & SUMMARY 

In this paper we have investigated whether spatially 
unresolved (integrated) spectroscopy of a spiral galaxy 
can be used to infer its mean interstellar oxygen abun- 
dance, as measured from observations of individual H 11 
regions. Our analysis makes use of a new integrated 
spectrophotome tric survey of nearby star -forming galax- 
ies presented bv iMoustakas fc KennicutH {2006), and ob- 
servations of more than 250 H il regions culled from the 
literature. Our final sample consists of 14, predominantly 



late-type spiral galaxies spanning a wide range of _B-band 
luminosity, inclination angle. Ha surface brightness, and 
dust content (Table For the first time, these data 
allow us to investigate the systematic effects of diffuse- 
ionized gas emission, dust reddening, and the slope of the 
abundance gradient on the determination of integrated 
abundances. 

We derive integrated oxygen abundances us- 
ing the emission-line meas urements presented by 
iMoustakas fc Kennicutt (2006), who carry out a careful 
treatment of underlying stella r abs orption. We apply 
both a theoretical (McGaug3 IT99lf) . and an empirical 
ijPilvugin fc ThuanI i2005.'l strong-line (i?23) abundance 
calibration to the observed emission-line fluxes, the 
reddening-corrected fluxes, and the emission-line equiv- 
alent widths (following Kob ulnickv fc P hillips 2003). 
We analyze the H ii-region data compiled from the 
literature self-consistently; in particular, we determine 
the H ii-region abundances using the same abundance 
calibration applied to the integrated emission lines, and 
we compute the de-projected galactocentric position 
of each object using robust inclination and position 
angles determined in the near-infrared. With these 
measurements, we derive the form of the abundance 
gradient for each galaxy in our sample (Fig. Table 121. 

Our principal result is that the integrated abun- 
dance of a normal disk galaxy correlates well with 
its characteristic gas-phase abundance, or the H 11- 
region abundance measured at p = 0.4 p25 (Figs. El and 
The typical scatter in the correlation is ±0.1 dex, 
independent of w h ether the jMcGaugh (199ij) or the 
iPilvugin fc ThuanI pOOl calibration is used (Table E)). 
We note, however, that the scatter in the abundance 
residuals is generally higher when using the Pilyugin fc 
Thuan calibration. These results confirm and extend 
the analysis carried out by Kobulnickv ct al. (1999) and 
IPilvugin et al.l l(2004al) . who used simulated integrated 
galaxy spectra. Integrated abundances based on the ob- 
served emis sion-line fluxes, or the emiss ion-line equiva- 
lent widths ijKobiilnickv fc Phillip^l2003l) are susceptible 
to additional systematic effects of order 0.05 — 0.1 dex 
with respect to the abundances based on the reddening- 
corrected fluxes, at least for the range of reddenings 
and stellar populations spanned by our sample. The 
systematic residuals using the equivalent-width method, 
however, can be minimized or eli minated by adopting a 
different value of a (see 21 and iKobulnickv fc Phillioij 
l2f)nl . 

We studied the residuals between the integrated and 
characteristic abundances of the galaxies in our sample 
to identify a physical origin for the source of the scat- 
ter. We found no signiflcant correlations between the 
residuals and the slope of the abundance gradient, the 
dust reddening, the inclination angle, or the Ha surface 
brightness (Figs. El and El and Table EJ. We note, how- 
ever, that Trcmonti ct al. (2004) found that inclination 
can signiflcantly bias integrated abundances, in the sense 
that edge-on galaxies may appear, on average, 0.2 dex 
more metal-poor than face-on galaxies. 

By comparing the emission-line sequences of galaxies, 
H II regions, and regions of diffuse-ionized gas, we found 
that the integrated [N 11] /Ha and [S 11] /Ha ratios of 
galaxies lie along the upper envelope of the sequence de- 
fined by high surface brightness H 11 regions, consistent 
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with a modest contribution from diffuse-ionized gas emis- 
sion to our integrated spectra (Fig. Cj). By comparison, 
the integrated [O iii]/H/3 and [O ii]/H/3 ratios of galaxies 
are relatively insensitive to contributions from diffuse- 
ionized gas emission, suggesting that the R23 parameter 
is a more robust abundance diagnostic for galaxies than 
strong- line calibrations that rely on either [S ll] or [N li] . 

Finally, we found that the reddening-corrected inte- 
grated abundances of the galaxies in our sample were 
equal to the gas-phase abundance at normalized galac- 
tocentric radii ranging from 0.2 to 0.8, depending on 
the slope of the abundance gradient. On average, the 
two abundances were equal at p/p25 — 0.47 ± 0.15, in- 
dicating that the commonly adopted radius used to de- 
fine the characteristic abundance of a disk galaxy (i.e., 
p/P25 — 0.4) is a reasonable assumption. 

We conclude that the integrated R23 parameter of star- 
forming disk galaxies is a robust tracer of the gas-phase 
oxygen abundance, even in the presence of a variety 
of systematic effects. The strength of the H/3, [O 11], 



and [O III] lines under typical nebular conditions, and 
their accessibility to ground-based optical spectrographs 
to z ^ 1, ensures that the R23 parameter will continue 
to be used to measure the build-up of heavy elements in 
the interstellar media of galaxies with cosmic time. 
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TABLE 1 
Integrated Oxygen Abundances 



Mb P25 i e log[E{H«)] E{B-V) 12+log(0/H)M9i" 12+log(O/H)pT05" 

Galaxy Type (mag) (arcmin) (deg) (deg) (erg s ^ pc ^) (niag) Observed Corrected EWs Observed Corrected EWs 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



NGC 1058 


SA(rs)c 


-18.4 


1.51 


26 


95 


32.6 


0.25 ± 0.08 


8.92 ± 0.03 


8.86 ± 0.03 


8.78 ±0.05 


8.34 ±0.04 


8.23 ± 0.04 


8.11 ±0.06 


NGC 2541 


SA(s)cd 


-18.2 


3.15 


64 


175 


32.0 


0.08 ± 0.08 


8.59 ± 0.04 


8.56 ± 0.04 


8.50 ±0.04 


8.28 ±0.04 


8.24 ±0.05 


8.15 ±0.06 


NGC 2903 


SB(s)d 


-20.2 


6.30 


59 


13 


32.3 


0.28 ± 0.08 


9.03 ± 0.02 


9.00 ±0.03 


8.95 ±0.04 


8.61 ±0.05 


8.52 ± 0.06 


8.42 ± 0.08 


NGC 3198 


SB(rs)c 


-19.9 


4.26 


72 


40 


32.1 


0.24 ±0.09 


8.74 ±0.06 


8.67 ±0.07 


8.68 ±0.06 


8.32 ±0.07 


8.19 ± 0.09 


8.24 ±0.08 


NGC 3344 


(R)SAB(r)bc 


-19.2 


3.54 


11 


150 


32.4 


0.15 ± 0.08 


8.85 ± 0.03 


8.82 ± 0.04 


8.72 ±0.05 


8.37 ±0.05 


8.30 ±0.05 


8.15 ±0.06 


NGC 3351 


SB(r)b 


-19.6 


3.71 


28 


163 


32.4 


0.55 ± 0.11 


> 9.05 


> 8.99 


> 8.98 


> 8.65 


> 8.48 


> 8.48 


NGC 3521 


SAB(rs)bc 


-20.1 


5.48 


61 


166 


32.6 


0.45 ± 0.07 


8.93 ± 0.02 


8.83 ± 0.04 


8.73 ±0.05 


8.48 ±0.03 


8.27 ±0.04 


8.13 ±0.05 


NGC 4254 


SA(s)c 


-20.3 


2.69 


18 


24 


33.2 


0.43 ± 0.07 


9.00 ±0.01 


8.94 ±0.02 


8.95 ±0.02 


8.53 ±0.02 


8.37 ±0.03 


8.42 ± 0.03 


NGC 4303 


SAB(rs)bc 


-20.0 


3.23 


42 


20 


32.9 


0.30 ±0.07 


8.99 ± 0.01 


8.95 ± 0.02 


8.96 ±0.02 


8.57 ±0.02 


8.46 ± 0.03 


8.50 ±0.03 


NGC 4321 


SAB(s)bc 


-21.0 


3.71 


39 


108 


32.7 


0.69 ± 0.09 


> 9.04 


> 8.99 


> 9.01 


> 8.77 


> 8.60 


> 8.67 


NGC 4651 


SA(rs)c 


-20.3 


1.99 


50 


80 


32.6 


0.24 ±0.07 


8.86 ± 0.02 


8.81 ± 0.03 


8.67 ±0.04 


8.45 ±0.03 


8.33 ± 0.04 


8.12 ±0.05 


NGC 4713 


SAB(rs)d 


-18.6 


1.35 


35 


50 


33.0 


0.13 ± 0.07 


8.71 ±0.02 


8.67 ±0.03 


8.65 ±0.03 


8.28 ±0.03 


8.21 ± 0.04 


8.20 ±0.04 


NGC 4736 


(R)SA(r)ab 


-19.4 


5.61 


38 


85 


32.2 


0.22 ± 0.07 


8.88 ±0.02 


8.83 ± 0.02 


8.63 ±0.03 


8.41 ±0.02 


8.30 ±0.03 


8.00 ± 0.04 


NGC 5194 


SA(s)bc 


-20.7 


5.61 


20 


58 


32.9 


0.39 ± 0.07 


9.01 ± 0.01 


8.96 ± 0.02 


8.95 ±0.02 


8.59 ±0.03 


8.45 ± 0.03 


8.45 ± 0.03 


Note. — Col. 


(1) Galaxy name; 


Col. (2) Morph. 


:>logioal 


tviDE fromlde 


Vaucoule, 


ara et alJITSSm: Col, 


(3) Absolute B 


-band magnitude. 


corrected for fc 


iroground Galactic c 


!xtinction (R^ = 


3.1: ISchlcffol ct al 


[^1996: O'Doniiell 



115541) . based on the apparent magnitude and distance gi ven by lMoustajtas^ T^ei ^ 12553) ; Col- (4) Radius of the major axis at the B25 mag arcsec" ^ isophote from l^TTaucouleurs et all 115511) ; Col. (5) Galaxy inclination 

angle based on t he JCs-band major-to-min or axis ratio rnT^^^^fB^POTn^TOcir Trdescribed in mTl\ except for NGC 5194, which is from lTTTTT^ |1{)74|) ; Col. (6) Galaxy position angle , measured positive from North to East, 
in t he Kg band ^^ ^ett et aH255(inro3) . as described iT" !^^ Col, f 71 Surface. Col. (8) Nebular reddening determined from the observed Ha/H/3 ratio IMoustaltas L kennicutj^ro^ . assuming an intrinsic ratio of 2.86 and 
the lO'l^onneli iW^W ^Alkv Way extinction curve; Oxygen abundances have been computed using the observed line fluxes [Cols. (9) and (12)], the reddening-corrected line fluxes [Cols. (10) and (13)], and the emission-line 
equivalent widths (EWs) [Cols. (11) and (14)]. 



^Oxygen abundances derived using either the lMc^aughl [T??! . M91) or the lPilvugln SU: TiTS^ HU^ , PT05) abundance calibrations. See iTTl for details. 
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TABLE 2 

H II-Region Oxygen Abundance Gradients 



12+log(U/H) 12+log(U/H) Gradient 12+log(U/H) 12+log(0/H) Gradient 

Galaxy at p = at p = 0.4p25 (dex P2k) at p = at p = 0.4p25 (dex p^r^) N(H ll) Refs. 

(1) (2) (3) (4) (5) (6) (7) ' (8) (9) 



NGC 1058 


9.09 ±0.01 


8.93 ± 0.01 


-0.41 ±0.02 


8.37 ±0.02 


8.35 ±0.01 


-0.05 ± 0.03 


6 


1 


NGC 2541 


8.46 ± 0.08 


8.51 ± 0.04 


0.13 ±0.17 


8.15 ± 0.19 


8.21 ±0.08 


0.16 ± 0.41 


19 


2 


NGC 2903 


9.11 ±0.01 


8.90 ±0.01 


-0.54 ±0.03 


8.56 ± 0.03 


8.37 ±0.01 


-0.46 ± 0.07 


45 


2,3,4,5 


NGC 3198 


8.98 ±0.07 


8.72 ± 0.03 


-0.65 ±0.13 


8.54 ±0.16 


8.32 ±0.07 


-0.54 ±0.30 


14 


2 


NGC 3344 


9.09 ±0.05 


8.76 ± 0.02 


-0.82 ±0.11 


8.34 ±0.10 


8.28 ±0.04 


-0.16 ± 0.22 


15 


2,4,6 


NGC 3351 


9.17 ±0.01 


9.09 ± 0.01 


-0.21 ±0.02 


8.68 ± 0.16 


8.60 ±0.07 


-0.21 ± 0.22 


14 


4,7,8 


NGC 3521 


9.17 ±0.06 


8.83 ± 0.03 


-0.84 ±0.20 


8.41 ± 0.14 


8.36 ±0.07 


-0.11 ± 0.50 


12 


2,7 


NGC 4254 


9.18 ±0.01 


8.99 ± 0.01 


-0.47 ±0.03 


8.62 ± 0.05 


8.46 ±0.02 


-0.40 ±0.07 


18 


4,9,10 


NGC 4303 


9.17 ±0.02 


8.87 ±0.01 


-0.76 ±0.05 


8.52 ± 0.05 


8.36 ±0.02 


-0.40 ±0.10 


22 


10,11 


NGC 4321 


9.26 ±0.03 


9.06 ± 0.01 


-0.48 ±0.07 


8.67 ±0.11 


8.50 ±0.04 


-0.42 ± 0.19 


9 


4,10 


NGC 4651 


8.94 ±0.01 


8.81 ± 0.01 


-0.32 ±0.03 


8.37 ±0.04 


8.32 ±0.02 


-0.13 ± 0.05 


7 


12 


NGC 4713 


9.03 ±0.05 


8.73 ± 0.01 


-0.76 ±0.10 


8.55 ± 0.08 


8.33 ±0.02 


-0.55 ± 0.17 


4 


12 


NGC 4736 


8.91 ±0.01 


8.85 ± 0.01 


-0.17 ±0.06 


8.44 ± 0.03 


8.31 ±0.04 


-0.31 ± 0.17 


15 


4,7,8 


NGC 5194 


9.20 ±0.01 


9.04 ± 0.00 


-0.41 ±0.03 


8.67 ±0.04 


8.54 ±0.01 


-0.33 ±0.08 


34 


4,7,13,14,15 



References. — tl) IPerEuson et all ('21 Igar'tslty et all Oj IBresolin et all f4) IMcCall et all 1155^1 : fSl lvan gee et all 11553 ); (Q) IVilcliez et all 

lliresoilii eT^fTBWHr rmi TCS el all lijiit J h (15l lUarneU all lUitij^ll . 

Note. — Col. (1) Galaxy name; Cols. (2) and (5) Central abundance (at radius p = 0) based on the derived abundance gradient; Cols. (3) and (6) Characteristic 
abundance (at p = 0.4p25), based on the derived abundance gradient; Cols. (4) and (7) Slope of the abundance gradient; Col. (8) Number of H II regions; Col. (9) 
H Il-region references. 

^Oxygen abundances and abundance gradients derived using either the |Mc(jausU |1{){)U . M91) or the lPilyugin lU: Tkuail llSOOgl PT05) abundance calibrations. See 11 
for details. 



TABLE 3 

Integrated vs. Characteristic Abundance Residuals 





M91'' 


PT05'' 


Method^' 


(Alog(0/H)> Alog(0/H)n,edian 


(Alog(0/H)) Al0g(0/H)^<3dian 


Observed 
Corrected 
EWs 


0.049 ± 0.055 0.047 
-0.003 ± 0.063 -0.008 
-0.063 ± 0.088 -0.042 


0.085 ± 0.083 0.093 
-0.028 ± 0.093 -0.009 
-0.111 ±0.128 -0.087 



''Integrated abundances based on the observed emission-line fluxes, the reddening- 
corrected line- fluxes, or the emission-line equivalent widths (EWs). See il3.1l for 
details. 

^Oxygen abundances d erived using either the IMcGauehl itTool M91) or the 
IPilvugin fc Thuaiil i2005l PT05) abundance calibrations. See >l3.1l for details. 



TABLE 4 

Residual Correlation Coefficients'' 



Variables 






M91' 




PT05'' 






Coetticient Probability 


Coetticient Pro 


bability 


Alog(0/H)obs vs. 


Slope 




-0.26 


0.42 


-0.07 


0.83 


Alog(0/H)„b, vs. 


E{B - 


V) 


0.08 


0.81 


0.31 


0.33 


Alog(0/H)obs vs. 


i 




0.42 


0.17 


0.22 


0.50 


Alog(0/H)obs vs. 


E{Ha) 




-0.39 


0.21 


-0.14 


0.66 


Alog(0/H)cor vs. 


yiope 




-0.22 


0.50 


0.16 


0.62 


Al0g(0/H)cor vs. 


E(B - 


V) 


-0.10 


0.76 


-0.13 


0.70 


Al0g(0/H)cor vs. 


i 




0.37 


0.24 


0.04 


0.90 


Al0g(0/H)cor VS. 


E(HQ!) 




-0.41 


0.18 


-0.20 


0.54 


Alog(0/H)Ew vs. 


yiope 




-0.29 


0.35 


-0.13 


0.68 


Alog(0/H)Ew vs. 


E{B - 


V) 


-0.02 


0.95 


0.19 


0.56 


Alog(0/H)Ew vs. 


i 




0.27 


0.39 


0.16 


0.62 


Alog(0/H)Ew vs. 


S(Ha) 




-0.11 


0.73 


0.15 


0.65 



''Spearman rank correlation coefficients and the probability that the variables 
are uncorrelated. We compare the residuals, Alog(0/H), for all three integrated 
abundance measurements (see HH.ll for details), against the slope of the abundance 
gradient, the dust reddening, E{B~V), the inclination angle, i, and the Ha surface 
brightness, E(HQf). 

^Oxygen abundances d erived using either the iMcGaueh (19911, M91) or the 
IPilvugin fc ThuanI ll2005l . PT05) abundance calibrations. See for details. 



